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Abstract 

We report on a numerical study of quantum transport in disordered two dimensional graphene 
and graphene nanoribbons. By using the Kubo and the Landauer approaches, transport length 
scales in the diffusive (mean free path, charge mobilities) and localized regimes (localization lengths) 
are computed, assuming a short range disorder (Anderson-type). In agreement with localization 
scaling theory, the electronic systems are found to undergo a conventional Anderson localization in 
the zero temperature limit. Localization lengths in weakly disordered ribbons are found to differ 
by two orders of magnitude depending on their edge symmetry, but always remain several orders 
of magnitude smaller than those computed for 2D graphene for the same disorder strength. This 
pinpoints the role of transport dimensionality and edge effects. 

PACS numbers: 73.63.-b,72.15.Rn,81.05.Uw 
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Recently, single graphene sheet could be isolated either from chemical exfoliation of bulk 
graphite fl. or by epitaxial growth on metal substrates through thermal decomposition 
of SiC [2|. These technological achievements have opened unprecedented opportunities to 
explore quantum transport in low dimensional carbon-based disordered systems ^, J]. 

Because of the unique electronic properties of the 2D graphene (massless Dirac fermions 
with linear dispersion and electron- hole symmetry), disorder effects and transport properties 
turn out to be unconventional. Theoretically, it has been shown that for long range impurity 
potentials, intervalley K ^ K' scattering between the two Dirac nodes could be strongly 
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reduced, resulting in anomalously low backscattering rates [5], extremely large elastic mean 
free paths and vanishingly small localization effects 6j. In contrast, for short range impurity 
potentials (where all types of scattering between K and K' are allowed), stronger quantum 
interferences could develop, leading to weak localization, or strong Anderson localization 
in the zero-temperature limit To date magnetotransport experiments either performed 
on exfoliated or epitaxial graphene have reported both weak antilocalization and weak lo- 
calization effects [Sj, confirming the sensitivity of 2D transport in graphene to the external 
random potential, whose precise origin remains unknown. 

Beyond 2D graphene physics, the transport properties of quasi- ID graphene nanorib- 
bons (GNRs) with width down to a few tens of nanometers have been characterized j^. 
In contrast to 2D graphene, the electronic properties of GNRs are strongly dependent on 
confinement effects and edge symmetries lOj. These new structures share similarities with 
carbon nanotubes, often viewed as rolled single graphene ribbons, and that have provided 
unique materials for investigating ID transport phenomena such as Luttinger liquid and 
Kondo physics or Anderson localization 

The issue of localization in graphene-based materials is currently highly debated from 
a theoretical standpoint. For instance, the measurement of a minimal conductivity of ~ 
2 — be^lh in samples for which charge mobilities change, however, by almost one order 
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of magnitude remains to be fully understood pj, 121, 13[. Indeed, a conventional treatment 
of disorder effects within the self consistent Born approximation (SCBA) yields o"™"^ ~ 
Ae'^/{hTx) for the two Dirac nodes {h the Planck constant), hence typically smaller by a 
l/vr factor with respect to the experimental data. Depending on the disorder model, the 



use of the Kubo approach suggest several scenarii to understand such discrepancies 



12|. 



Besides, the role played by quantum interferences and the transition to a localization regime 
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in graphene and GNRs remain poorly explored but fiercely debated 

In this Letter, by using both the Kubo and Landauer approaches, the transport length 
scales in 2D graphene are investigated and compared with those of the quasi- ID GNRs. The 
disorder (Anderson-type) is introduced via random fluctuations of the onsite energies of the 
vr-orbitals, which mimics a short range scattering potential that has been widely studied in 



the past as a generic disorder model in the framework of localization theory |13l. Il4l|. For 2D 
graphene, a real space order Kubo method [l6i] is used to compute the energy-dependent 



elastic mean free path (fg), charge mobilities (/x) and semiclassical conductivities (cTsc) in 
the diffusive regime, before quantum interferences come into play. Beyond the diffusive 
regime, the energy-dependent localization length (^) is extracted from the analysis of the 
transition from weak to strong localization, following the scaling theory phenomenology . 
Quantum transport in GNRs with different chiralities (ziz-zag and armchair type) and same 
disorder potential is also investigated within a Landauer approach For GNRs widths 
in the range ~ 20 — 80 nm (within the experimental scope 9|), it is found that edge effects 
strongly enhance the impact of disorder, which results in localization lengths several orders 
of magnitude smaller than those obtained in 2D graphene for the same disorder strength. 
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FIG. 1: (color online) (a): DoS of an ideal (dashed lines) and for disordered graphene sheets for 
several values of W = 1, 1.5, 2, 2.5 (b): Zoom in the energy area around the CNP. 



The low energy electronic properties of 2D graphene are accurately described by the vr- 
orbital tight-binding hamiltonian, which is a first nearest neighbor two centers orthogonal Pz 
model, with onsite energies Ec = eV for all orbitals and hopping term 70 = 2.7 eV. To mimic 
short range disorder, a white noise uncorrelated Anderson type disorder is introduced as a 
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random fluctuation of tlie onsite energies of the hamiltonian {e = ec + Se). The scattering 
potential can thus be characterized by a single parameter W which defines the range of 
energy variations {6e G [—W'yo/2, +W'yo/2]), and thus allows to tune the disorder strength. 
In what follows W = [0.5, 2.5] enables the exploration of all transport regimes taking place 
in disordered 2D graphene and GNRs. 

In Fig. [H the Density of States (DoS), computed with a Lanczos-type method [l6|, is 
reported as a function of disorder strength. The disorder-free DoS (dashed line) shows the 
typical behavior with a linear increase at low energy and the presence of two sharp Van 
Hove singularities at E = ±70. As W is increased, two opposite behaviors are observed. At 
high energies. Van Hove singularities are smoothen whereas close to the charge neutrality 
point (CNP), disorder enhances the DoS in agreement with prior analytical results [51] (see 
Fig. [T]d for a close-up). 

To investigate quantum transport in the 2D disordered graphene, an efficient real space 
and order N Kubo method is employed 16|. In this formalism, the mean free path ie{E), 
the semiclassical conductivity asc{E) and the charge carrier mobility fi{E) are deduced 
from the energy and time dependence of the diffusion coefficient D{E,t) = {AR'^){E,t)/t 
(where {AR'^){E,t) is the quadratic spread of random phase wavepackets propagated in the 
graphene sheet). 

In Fig. [2]3, the time dependence of D{E,t) at the CNP and at E = 0.1 eV are reported 
for two values of W ITj. Different transport regimes follow each other as a function of the 
propagation time (or length). As expected, D{E,t) first scales linearly with t at short times 
owing to the absence of elastic scattering. This linear scaling is followed by a saturation 
of -D at a maximum value D^ii^{E), that pinpoints the occurrence of a diffusive regime for 
which D ~ vie (with v a group velocity and ie the elastic mean free path [16]). As evidenced 
in Fig. the saturation time decreases with increasing disorder strength (W) or increasing 
charge energy (E). At longer times, D{E,t) decreases owing to quantum interferences effects 
and localization phenomena [l^ . The full energy-dependence of ie is given in Fig. [2b for 
increasing W. 

The strong enhancement of ie around the CNP results from the cusp in the DoS, which 
implies a reduced number of scattering processes. However, ie{E) drops from 180 nm at 
W = 1.0 to 10 nm a.t W = 2.0, as a consequence of the increase of the DoS close to CNP. 
In the weak disorder case {W G [0.2 — 0.7] not shown here), when the DoS at the CNP is 
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FIG. 2: (color online) (a): Energy-dependent charge mobility for several values of W. (b): ie for the 
same disorder strengths. The legend is the same as in Fig. [TJ (c): Diffusion coefficient D{E,t) as 
a function of time for various disorder strengths and Fermi energies. D{E,t) has been normalized 
with respect to its maximum value Dmax{E) to allow an easier comparison between the different 



curves. 



almost unchanged with respect to disorder- free graphene case, the behavior of ie{E) as a 
function of W is in good agreement with the Fermi golden rule (FGR) i.e. ie{E) oc l/W"^. 
For higher values of W (1.0 to 2.0) slight deviations to the FGR are expected, since the 
weak disorder approximation is not strictly applicable anymore. 

Fig. [2^ shows the corresponding charge mobilities deduced from fi{E) = asc{E) / en{E) , 
where cXsc = e'^p{E)v{E)i(. is the semiclassical conductivity deduced from the Einstein for- 
mula, p{E) is the DoS, n{E) is the charge density at energy E, and is e the elementary 
charge. The energy-dependence of yu(-E) and ie{E) are similar, and the sharp increase of 
uiE) in the vicinity of CNF is in good qualitative agreement with experimental observations 

Experiments show that the conductivity (down to a few Kelvin) is almost constant close 
to the CNP, crfi^^ = 0) ~ 3 — 5e^//i, and weakly dependent on the value of the charge 
mobility y, l^l. On the theoretical side, within the SCBA the semiclassical part of the 
conductivity due to short range disorder is found to be asc = 4e^/(/i7r) [5]. By using the 
Kubo formalism, it was further found that a{E = 0) strongly depends on the nature of the 
scattering potential (short or long range) [l^. Our numerical results are shown in Fig. ^jp. 
asc clearly remains larger or equal to amin = '^Gq/ti {Gq = 2e^//i), a fact that would be 
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FIG. 3: (color online) (a): (,{E) for three disorder strengths (same legend as in Fig. [T]). (b): 
Energy-dependent semiclassical conductivity for the same disorder strengths. 

consistent with Mott argument [l^ , although localization effects are further observed at all 
energies (see hereafter). Besides, the shape of the energy-dependence of conductivity is in 
perfect agreement with prior analytical results derived for short range disorder within the 
SCBA fl. 

However, the conductivity would not be sensitive to localization effects only in the pres- 
ence of some decoherence mechanisms such as electron-electron scattering of electron-phonon 
coupling In contrast, as previously seen in the time-dependence of the diffusion coef- 
ficient (Fig. [Jt), our zero-temperature calculations evidence the contribution of localization 
effects that dev elop beyond the diffusive regime. The 2D localization length ^ can be evalu- 
ated as follows [14'|: Whatever the disorder model, the quantum correction to the conductiv- 
ity is expected to scale as Acr(L) = (Go/n) \ia.{L/ie), where L is the length scale associated 
with the propagation time. The localization length ^ is given by Aa{L = C.) = cTsc, i-e. 
^ = £eexp(7r(Tsc/G'o). Our results are reported in Fig. [3k for several disorder strengths. The 
energy-dependence of ^ is mainly dominated by that of cr^c- As a result, although is 
strongly increasing as the Fermi level moves towards the CNP (undoped case), the behavior 
of ^ shows an opposite trend, with a minimum value at CNP. 

Recently, the possibility to fabricate quasi- ID graphene nanoribbons has opened new 
perspectives for future carbon-based nanoelectronics 2|, |9l]. It is thus important to evaluate 
the effects of disorder in this situation of lower dimensionality. Indeed, the bandstructures 
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FIG. 4: (color online) (a): Conductance for ideal ziz-zag (solid blue line) and armchair (dashed 
red line) GNRs with width of ~ 20 nm. (b): Conductance for a single disorder configuration of 
a zig-zag (solid blue line) and an armchair (dashed red line) CNR with width ~ 20 nm and for 
W = 0.5. (c): Configuration averaged (over ~ 400 samples) normalized conductance as a function 
of GNR length for both zig-zag and armchair GNRs. The solid blue (dashed red) arrow shows the 
energy at which the calculations for the zGNR (aGNR) have been performed, (d) and (e): Same 
informations as for (b) and (c) but for a larger disorder strength (W = 2). 



of ideal GNRs with width below ~ 100 nm and well defined edge symmetries (zig-zag or 
armchair types) are dominated by confinement effects and Van Hove singularities 10|] , simi- 
larly to carbon nanotubes Zig-zag type GNRs show very peculiar electronic properties 
with wavefunctions sharply localized along the GNRs edges at low energies, which will sig- 
nificantly influence their transport properties. In contrast, armchair-type GNRs share more 
similarities with metallic nanotubes. Some prior studies have already addressed the impact 



of homogeneous disorder on quantum transport in GNRs 181. 

nn 

By using a Landauer approach [15|, Il9|, ie{E) and ^{E) can be estimated from the scaling 
analysis of the average conductance of the GNRs jl9]. In particular, ^{E) can be accurately 
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extracted from the exponential decay of the conductance (averaged over ~ 400 configura- 
tions) versus length. 

In Fig. m the conductance for both types of symmetries are shown in the ideal case and 
for a single disorder configuration for weak {W = 0.5) as well as strong disorder {W = 2). 
The length scaling behaviors of the conductance and the localization lengths (for both types 
of GNRs edge symmetries) are respectively reported in Fig. HI: and and in Fig. [31 For 
disorder as large as W = 2, the localization lengths are similar for both types of ribbons, 
showing that edge symmetry does not play any role. In contrast, ^ is up to two orders 
of magnitude larger in zig-zag than in armchair GNRs in the low disorder limit {W = 
0.5), for all values of energy in the plateau around CNP (some illustrative ^ values at 
chosen energies are shown in Fig. [3]). This can be simply understood as the consequence 
of the lower dimensionality of transport in the case of zig-zag symmetry, driven by more 
confined wavefunctions [10]. However, ^ remains always several orders of magnitude smaller 
in GNRs than in 2D graphene, whatever the disorder strength. Similar results are obtained 
for GNRs with larger width of ~ 80 nm (not shown). This points out the predominant role 
of dimensionality and edge effects on quantum transport especially in the low disorder limit. 

In conclusion, by studying the transport properties in both disordered 2D graphene and 
GNRs (for short range scattering potential), the impact of edge symmetries and transport 
dimensionality on the transport length scales was outlined. Despite the simplicity of the 
Anderson model, some of the reported transport features may be generic to other types of 
disorder such as chemical doping, surface functionalization, or topological defects. However, 
the contribution of long range potentials (e.g. due to ionized impurities trapped in an oxide) 
deserves further consideration. 
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